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ABSTRACT 

We study, in this paper, the non-Gaussian features of the mass density field 
of neutral hydrogen fluid and the Lya transmitted flux of QSO absorption spec- 
trum from the point-of-view of self-similar log-Poisson hierarchy. It has been 
shown recently that, in the scale range from the onset of nonlinear evolution to 
dissipation, the velocity and mass density fields of cosmic baryon fluid are ex- 
tremely well described by the She-Leveque's scaling formula, which is due to the 
log-Poisson hierarchical cascade. Since the mass density ratio between ionized 
hydrogen to total hydrogen is not uniform in space, the mass density field of 
neutral hydrogen component is not given by a similar mapping of total baryon 
fluid. Nevertheless, we show, with hydrodynamic simulation samples of the con- 
cordance ACDM universe, that the mass density field of neutral hydrogen, is also 
well described by the log-Poisson hierarchy. 

We then investigate the field of Lya transmitted flux of QSO absorption 
spectrum. Due to redshift distortion, Lja transmitted flux fluctuations are no 
longer to show all features of the log-Poisson hierarchy. However, some non- 
Gaussian features predicted by the log-Poisson hierarchy are not affected by the 
redshift distortion. We test these predictions with the high resolution and high 
S/N data of quasars Lya absorption spectra. All results given by real data, 
including /5-hierarchy, high order moments and scale-scale correlation, are found 
to be well consistent with the log-Poisson hierarchy. We compare the log-Poisson 
hierarchy with the popular log-normal model of the Lya transmitted flux. The 
later is found to yield too strong non-Gaussianity at high orders, while the log- 
Poisson hierarchy is in agreement with observed data. 
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1. Introduction. 

Baryon matter of the universe is mainly in the form of intergalactic medium (IGM), of 
which the dynamics can be described as compressible fluid. Luminous objects are formed 
from baryon matter in the gravitational well of dark matter. Therefore, the dynamical state 
of baryon fluid in nonlinear regime is crucial to understand the formation and evolution of 
the large scale structures of the universe. In the linear regime, the baryon fluid follows the 
mass and velocity fields of coUisionless dark matter. In the nonlinear regime, however, the 
baryon fluid statistically decouples from the underlying dark matter field. The statistical 
behavior of baryon fiuid is no longer described by a similar mapping of the underlying dark 
matter field (e.g. Pando et al. 2004). 

It was first pointed out by Shandarin and Zeldovich (1989) that the dynamical behavior 
of baryon matter clustering on large scales is similar to turbulence. The expansion of the 
universe eliminates the gravity of uniformly distributed dark matter. The motion of baryon 
matter on scales larger than dissipation is like that of matter moving by inertia. In this 
regime, the evolution of baryon matter is scale-free and dynamically like fully developed 
turbulence in inertial range. The turbulence of incompressible fluid leads the energy passes 
from large to the smallest eddies, and finally dissipates into thermal motion. While the 
clustering of cosmic baryon fiuid is also due to the transform of density perturbations on 
different scales, and finally falls and dissipates into virialized halos of dark matter. Yet, the 
turbulence of incompressible fiuid is rotational (Landau & Lifshitz 1987), while the clustering 
of cosmic matter is irrotational, because vorticities do not grow in an expanding universe 
(Peebles 1980). 

Nevertheless, the turbulence-like behavior of cosmic baryon fiuid has been gradually 
noticed. First, the dynamics of growth modes of the cosmic matter is found to be sketched 
by a stochastic force driven by Burger's equation (Gurbatov, et al 1989, Berera & Fang, 
1994). The Burger's equation driven by the random force of the gravity of dark matter 
can also sketch the evolution of baryon fluid, if cooling and heating are ignored (Jones 
1999; Matarrese & Mohayaee 2002). Later, the Burger's fluid is found to show turbulence 
behavior if the Reynolds number is large enough (Polyakov 1995; Lassig 2000; Bee & Frisch 
2000; Davoudi et al. 2001). The Reynolds number of IGM at nonhnear regime actually is 
large. Therefore, we may expect that, in the scale free range, the dynamical state of cosmic 
baryon fluid should be Burger's turbulence. The turbulence of Burger's fluid is different from 
the turbulence of incompressible fluid. The later consists of vortices on various scales, while 
the former is a collection of shocks. 

With the cosmological hydrodynamic simulation based on Navier-Stokes equations in 
which heating and cooling processes are properly accounted, it has been found that the 
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velocity field of the IGM consists of an ensemble of shocks, and satisfies some scaling relations 
predicted by Burger's turbulence (Kim ct al. 2005). This result reveals that the turbulence 
features of cosmic baryon fluid are independent of the details of dissipation (heating and 
cooling) mechanism if we consider only the scale free range, i.e. from the scale of the onset 
of nonlinear evolution to the scale of dissipation, say Jeans length. 

A new progress is to show that the velocity field of cosmic baryon fluid can extremely 
well described by She-Leveque's (SL) scaling formula (He ct al. 2006). The SL formula is 
considered to be the basic statistical features of the self-similar evolution of fully developed 
turbulence. Very recently, the non-Gaussianitics of mass density field of the hydrodynamic 
simulation samples are found to be well consistent with the predictions of the log-Poisson 
hierarchy, which originates from some hidden symmetry of the Navier-Stokes equations. 
This hierarchical model gives a unified explanation of non-Gaussian features of baryon fluid, 
including the intermittence, hierarchical relation, scale-scale correlations etc (Liu & Fang 
2008). These results strongly indicate that, in the scale free range, dynamical state of 
cosmic baryon fluid is similar to a fully developed turbulence. 

In this paper, we investigate the log-Poisson hierarchy of cosmic baryon fluid with 
observed data - the Lya transmitted flux of quasar absorption spectrum, which is due to 

the absorptions of quasar continuum by the diffusely distributed neutral hydrogen (Bi et al 
1995, Bi & Devidscn 1997, Ranch 1998). These samples offer a unique way to study the non- 
Gaussian feature of cosmic baryon ffuid. It has been known for a long time that the fields 
of Lja forests and transmitted fiux arc highly non- Gaussian. Observation samples of Lya 
forests and transmitted fiux show scale-scale correlation (Pando et al 1998), intermittence 
(Jamkhedkar et al. 2000; Pando et al. 2002; Feng et al. 2003), non-thermal broadening of 
H I and He II Lya absorption fines (Zheng et al. 2004; Liu et al. 2006). We will show that 
the log-Poisson hierarchy provides a crux to understand the non-Gaussian behavior. 

The outline of this paper is as follows. §2 gives an introduction of the log-Poisson 
hierarchy. §3 shows that the neutral hydrogen component of cosmic baryon fluid is of log- 
Poisson hierarchy. In §4, we study the log-Poisson hierarchy of the held of Lya transmitted 
flux with observed samples of quasar absorption spectra. A comparison between log-Poisson 
hierarchy and log-normal model is also presented in §4. The conclusion and discussion are 
given in §5. 
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2. Log-Poisson hierarchy 



2.1. 



Structure function 



To describe the statistical properties of an isotropic and homogenous random field p(x) , 
it generally uses correlation functions of 5p(x) = p(x) — p, p being the mean of density 
For instance, a two-point correlation function is ((5p(x)(5p(x')). To reveal the turbulence-fike 

behavior, we use variable 6pr = p(x + r) — p(x), where r = |r|. The variable Spr is very 
different from variable 5p(x). The later can be larger than p, but cannot be less than —p, 
and therefore, for a nonlinear field, the distribution of 5p(x) generally is skew; however, 
the distribution of 6pr is symmetric with respect to positive and negative 6pr if the field is 
statistically uniform. 

With 5p(x) the clustering and non-Gaussianity of mass density p(x) are measured by 
two and multiple point correlation functions of Sp{x). With Spr, the statistical features are 
described by structure function defined as 



where p is the order of statistics, and the average (...) is taken over the ensemble of density 
field. A comparison of the correlation function and structure function has been analyzed in 
detail by Monin & Yaglom (1975). The 2nd structure function 5*2 = (|(^PrP) as a function r 
(scale) actually is the power spectrum of the mass density field p(r) (Fang & Feng 2000). 

In the scale-free range of the dynamical equations and initial conditions, the structure 
function is scale-invariant, and therefore, it is generally expressed as a power law of r 



^{p) is called intermittent exponent. Since the pioneer work of Kolmogorov (1941), it is 
believed that the relation of ^(p) vs. p is related to the scale-covariance of the dynamical 
equations and initial conditions. For fully developed turbulence of Navier-Stokes fluid, ^ (p) 
is a nonlinear function of p. Since then many hierarchy models for interpreting ^ (p) have 
been proposed (Frisch 1995). Finally the best model is given by the SL scafing formula (She 
& Leveque 1994), which is yielded from the Log-Poisson hierarchy process (DubruUe 1994). 

Although the cosmic baryon fluid is not incompressible, samples of mass and velocity 
fields of cosmic baryon fluid produced by the cosmological hydroidynamic simulation of the 
concordance ACDM model show in good arrgement with the SL scafing and log-Poisson 
hierarchy. This is not surprising, because the hierarchical structure model is mainly based 
on the invariance and symmetry of nonlinear dynamical systems. Therefore, systems other 



S,{r) = {\Sprn, 



(1) 



Sp{r) oc r' 



(2) 
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than the Navier-Stokes incompressible fluid wiU also show the SZ scaling and log-Poisson 
hierarchy if they have the similar invariance and symmetry (She 1997). 

2.2. Log-Poisson hierarchical cascade 

The scenario of hierarchical clustering has been widely used to describe nonlinear evolu- 
tion of the mass fleld of cosmic matter. We will flrst give the basic assumptions of log-Poisson 
hierarchy cascade, and then discuss the physics behind this model. 

The log-Poisson hierarchy assumes that, in the scale free range, the variables (density 
fluctuation) \5pr\ on different scales r are related by a statistical relation as (DubruUe 1994, 
She & Waymire 1995), 

\5pr^\ = Wrir-2l<^Pnl> (3) 

where 

Wr,r, = r{ri/r2V. (4) 
In eqs.(3) and (4), ri > r2, and m is a random variable with the Poisson PDF as 

P{m)^eM-K,r,)K,rJm\. (5) 

To insure the normalization (H^rir2) = 1) where (...) is over m, the mean \r^r2 of the Poisson 
distribution should be 

A,,,, = 7[ln(n/r2)]/(l-/3). (6) 

The model eq.(3) describes how a density fluctuation \5pr^ \ on larger scale ri statistically 
related to fluctuation |5pr2l on smaller scale r2- The log-Poisson model depends only on the 
ratio r\jr2- Thus, it is scale invariant. The random variable m can be considered as the 
steps of the evolution from I^PrJ to |5pr2l- When /5 is smaller, only the evolutionary process 
with lower steps is important. 

For a Gaussian field, variables bpr^ and bpr^ have to be statistically independent. It likes 
that the Fourier modes with different wavenumber k\ oc 1/ri and k2 oc l/r2 are statistically 
independent. Therefore, a Gaussian field has to be /3 = 1. Thus, the parameters /3 < 1 is a 
measure of the deviation from Gaussian field. The meaning of 7 will be given later. 

Among hierarchical clustering models, the log-Poisson hierarchy has the following fea- 
tures. First, the relation between l^pr-J and |5pr2l given by eq.(3) is a multiplicative random 
process. A random multiplicative cascade generally yields a non-Gaussian field. That is, 
even the field on large scale ri is Gaussian, it will be non-Gaussian on small scale r^- This 
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is different from additive random process (e.g. Cole & Kaiser 1988), which generally yields 
Gaussian field (Pando et al. 1998). 

Second, the cascade from scale ri to r2, and then to r3 is identical to the cascade 
from ri to r^. It is because W^ir-s — ^r^r^Wrirz — l3^{i'i/fzV, where N is again a Poisson 
random variable with Xr^rs — ^^1^2 + Xr^rs- Therefore, the log- Poisson hierarchy removes an 
arbitrariness in defining the steps of cascade from ri to or r2 to r^- This arbitrariness is a 
major shortcoming of many hierarchical clustering models, one of them, for instance, is the 
clustering hierarchy models (Soneira & Peebles 1977, Peebles 1980). 

Third, although the log-Poisson hierarchical process is discrete in terms of the discrete 
random number m, it is infinitely divisible. That is. there is no lower limit on the difference 
Ti — r2- It can be infinitesimal. This is consistent with the continuous variable r used 
in the hydrodynamic equation of cosmic baryon fiuid. The infinite divisibility can not be 
modeled with hierarchy of discrete objects with finite size. Log-normal model is also infinitely 
divisible. However, their asymptotic behavior of ^(p) at larger p is unbound. We will compare 
the log-normal model with log-Poisson in §4.4. 



2.3. Intermitted exponent 

With log-Poisson hierarchy eq.(3), the intermittent exponent ^(p) is found to be (Liu & 
Fang 2008) 

e(p) = -7[p-(l-/5^)/(l-/3)]. (7) 

When /3 ^ 1, we have ^{p) = 0. Therefore, = 1 is a Gaussian field. A field with f3 < 1 
is called intermittent. Eq.(7) requires ^(1) = 0. Therefore, the difference of an intermittent 
field from Gaussian is mainly given by the term 7(1 — f3^)/ (1 — f3) of cq.(7). 

From eq.(7), the power spectrum 5*2 (r) = const is flat. This is not generally applicable 
for the cosmic density field. In the scale free range, the power spectrum of mass field is of 
power-law. Thus, we should generalize the log-Poisson hierarchy eq. (3) by replacing ^p^i and 
5pr2 with rfSpr^ and r25pr2- In this case, the intermitted exponent ^{p) is 

e(p) = -«P - 7[P - (1 - /9^)/(l - /9)]- (8) 

Eqs.(2) and (7) yield 5*2 (r) oc r~^", and therefore, the parameter 2a is the index of power 
spectrum. When a = 0, eq.(8) is the same as eq.(7). 
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2.4. /3-hierarchy 

Since \6pr\^ is the p-th moment of the 5pr, for high p, one can attribute the structure 
function Sp{r) to the events focated at the tail of the probabihty distribution function (PDF) 
oi 5pr. To pick up the structures, which dominant the p-order statistics of dpr, we define 

F,(r)^ Vi(r)/5,(r). (9) 

Prom equations (2) and (8), equation (9) gives 

Fp(r) = Ar-"-T(i-^"). (10) 

where the constant A is independent of r and p. Por an intermittent field /? < 1, we have 
Foo{f) = Ar~°'~'^. Thus, from eq.(lO), one can find 

Equation (11) is invariant with respect to a translation in p. Since Fp{r) measures the 
structures dominating the p order statistics, the larger the p, the larger the contribution 
of strong-clustered structures to Fp{r). Therefore, equation (11) describes the hierarchical 
relation between the stronger (or high p) and weaker (or low p) clustering. In the scale-free 
range where Fp+i{r) / F^{r) < 1, we have Fp{r)/F^{r) < Fp+i{r) / F^{r) if /3 < 1. That 
is, for an intermittent field, weak clustering structures are strongly suppressed with respect 
to the strong clustering; the smaller the (3, the stronger the suppression of weak clustering 
structures. 

Prom eq.(ll), we have 

lnFp+i(r)/F3(r) = /3\nFp{r)/F2{r). (12) 

This, for allr and p, \n[Fp^i{r) / F^lr)] vs. \n[Fp{r) / F2{r)] should be on a straight line with 
slope [3. It is called /9-hierarchy. Eq.(12) does not contain parameter 7 and term F^[r), and 
therefore, it is a testable prediction of the log-Poisson hierarchy. 

Therefore, log-Poisson hierarchy links the sizes of fluctuation structures [eq.(3)] as well 
as their amplitude or intensity [eq.(ll)]. This is different from hierarchy models, which give 
only the relationship between the sizes of objects. The hierarchical relation eq.(ll) actually 
is the origin of the log-Poisson hierarchy. That is, the hierarchical relation like eq.(ll) is the 
first recognized to be held for dynamical system described by the Naiver-Stokes equations, or 
equations close to the Naiver-Stokes fiuid (DubruUe 1994, Leveque & She 1997). Therefore, 
the SL formula has also been successfully applied to describe the mass fields of compressible 
fiuid (Boldyrev et al 2002, Padoan et al 2003). 



Fp+i{r) 
Fooir) 
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2.5. 7"''6l^t6d non-Gaussianities 

The ratio between higher order to second order moments {6pfF)/ (Spf.)'^ is a popular tool 
to measure non-Gaussianity. When p = 2, the ratio is kurtosis. For a Gaussian field, the 
ratio is independent of r, and equal to 

(13) 

For log-Poisson hierarchy, one can show(Liu & Fang 2008) 

In = Kp\nr + const ( 14) 

The moment ratio la{{6p^P)/ (Sp^y) is hnearly dependent on Inr (scale free) with the slope 
Kp given by 

As expected, for Gaussian field {j3 ^ 1), Kp = 0, i.e. the ratio of moments is independent 
of r [Eq.(14)]. That is, Kp only depends on the last term on the right hand side of eq.(8), 
regardless the parameter a. 



Other useful non-Gaussian detector is the scale-scale correlation defined as 

{Sp^.,SpP,) 

m){sp?^..y 

which describes the correlation between the density fiuctuations on scales ri and r2. Ob- 
viously, for a Gaussian field, C^^f^^ = 1. The clustering of cosmic large scale structure in 
the nonlinear regime essentially is due to the interaction between the modes of fluctuations 
on different scales (e.g., Peebles 1980). Therefore, it is important to detect the scale-scale 
correlation of cosmic baryon matter. 

If the ratio r2/ri is fixed, the log-Poisson hierarchy predicts the scale-scale correlation 
to be (Liu & Fang 2008) 

Cf;^,^ = Bir2/n)rfP^~'^^P\ (17) 
where the factor B{r2/ri) is constant when the ratio r2/ri is fixed. This is because the 
log-Poisson model is invariant with respect to scale dilation. From eqs.(7) or (8), we have 

e(2p) - 2^{p) = -7(1 - ppf/il - (3) (18) 

Thus, if r2/ri remains constant, the relationship of InC^^'^^^ vs. Inri should be a straight 
fine with slope —7(1 — /3^)^/(l — P). This slop is also only dependent on the last term on 
the right hand side of eq.(8), regardless the parameter a. 
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3. Log-Poisson hiercirchy of neutral hydrogen density field 

Lja absorption spectrum depends on the distribution of neutral hydrogen. We first 
study the mass density field of diffused neutral hydrogen. Since the UV background radiation 
is uniform and does not introduce special spatial scale, one may expect that the mass density 
field of neutral hydrogen is also a field of the log-Poisson hierarchy. However, the ratio of the 
neutral hydrogen density Phi(x) to total hydrogen density Ph(x) is not spatially constant, 
because the temperature-density relation is of multiphase (He et al 2004, 2005). Therefore, 
the mass field of the neutral hydrogen would not have the same log-Poisson hierarchical 
features as of p(x). 

3.1. Simulation samples 

The simulation samples of the fields of baryon fiuid are generated by the same way as 
Liu & Fang (2008), which is based on the hybrid hydrodynamic N-body code of Weighted 
Essentially Non-oscillatory (WENO) scheme (Feng et al. 2004). It is in Eulerian scheme 
and suitable to analyze the fluid in high as well as in low mass density areas. The WENO 
samples have been successfully applied to reveal the self-similar hierarchy behavior of cosmic 
baryon fluid (Kim et al. 2005; He et al. 2006; Liu & Fang 2008), to explain the HI and Hell 
Lyct absorption in quasar spectra (Liu et al. 2006), and to study the relation between X-ray 
luminosity and temperature of groups of galaxies (Zhang et al. 2006). It is also used for 
studying Lya leaks at high redshift (Liu et al. 2007). 

The simulation is performed in a comoving cubic box of 100 Mpc with 512'^ grids 
and an equal number of dark matter particles. The grid size is 100/512 ~ 0.2 h^^ Mpc, 
which is less than Jeans length, and therefore, it is enough to catch the statistical behavior 
in the scale-free range, i.e. larger than the Jeans length and less than non-linear scale. 
We use the concordance ACDM cosmology model with parameters Qni=0.27, Qb=0.044, 
nA=0.73, h—0.71, (J8=0.84, and spectral index n = 1. The transfer function is calculated 
using CMBFAST (Seljak & Zaldarriaga 1996). We take a primordial composition of H and 
He (X=0.76, F=0.24). 

The ionization fraction is calculated with ionization-recombination equilibrium under 

a uniform UV radiative background, of which the intensity is adjusted to flt the mean of 
observed Lja tarnsmitted flux. We produce the distributions of the mass density of hydrogen 
/?(x), the fraction of neutral hydrogen /hi(x), temperature T(x), and velocity v(x) at redshift 
z ~ 2.5, which is the mean redshift of observed samples used in §4. 

We randomly sample 10,000 one-dimensional sub-samples to simulate the Lya trans- 
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mitted flux (§4). To estimate the errors, we divided the 10,000 samples into 10 subsamples, 
each of which has 1,000 hne samples. 
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Fig. 1. — Function Fp(r) vs. r of simulation samples of the mass density fleld of neutral 
hydrogen in the scale range 0.9 < r < 15 h~^ Mpc, and p — 0.5 x n, with n — 1, 2. ..8 from 
bottom to top. 



3.2. Log-Poisson hierarchical statistics 

The variable Spr = p(x + r) — p(x) of density field p(x + r) can be calculated by a 
discrete wavelet transform (DWT) as 

Spr,i = J pm{x)ipj,i{x)dx, (19) 

where ipj^i{x) is the base of discrete wavelet transform (e.g. Fang & Thews 1998). For a 
one-dimensional sample in physical space from a; = to L, the scale index j is related to 
the scale r hy r = and the position index / is for the cell located at a; = to 

(/ + 1)L/2K Because the DWT bases iJj,i{x) are orthogonal, the variables Sp^^i do not cause 
false correlation. They are effective to describe turbulence of fiuid (Farge 1992). For a given 
scale r or j, the statistical average (...) of eq.(5) is over all cells I. We will use the Harr wavelet 
to do the calculation below. We also repeat the calculations with wavelet Daubechies 4. 

Figures 1-3 demonstrate that the density field Phi(x) of neutral hydrogen shows the 
features of log-Poisson hierarchy. Figure 1 shows the functions -Fp(r) of the field pni in 
the physical length scale range 0.9 < r < 15 h~^ Mpc and orders of p = 0.5 x n with 
n — 1,2. ..8. For all p, InFp(r) can be well fitted by a straight line of Inr. It is consistent 
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Fig. 2. — ln[Fp_|_i(r)/F3(r)] vs. ln[Fp(r)/F2(r)] for simulation samples of the mass density 
field of neutral hydrogen. The data points are on scales r =1.78, 3.57, 7.14 and 14.3 h~^ 
Mpc and the statistical order to be p = 1, 1.5, 2, and 2.5. 



with eq.(lO). Whenp > 3, -Fp(r) given in Figure 1 is almost independent of p. This indicates 
that Fp{r) for higher p. Thus, from eq.(lO), (3 should be less than 1, and therefore, 

the field is intermittent. In this case, the slope of the straight lines with higher p has to be 
equal to —a — 7. Figure 1 shows a + 7 = 1.0. 

Figure 2 is the /3- hierarchy eq.(12) of the density field pni, in which the error bars are 
given by the ranges of ln[Fp+i(r)/F3(r)] and ln[Fp(r)/F2(r)] of the 10 subsamplcs, each of 
which consists of 1,000 lines. The /3-hierarchy is held for all r in 1 < r < 15 h~^ Mpc and 
p =0.5 to 4. It yields (5 = 0.73 ± 0.19. 

Figure 3 plots the intermittent exponent ^(p) [eq.(8)]. It shows that ^(p) can be well 
fitted with eq.(8) with parameters [3 = 0.60, a + 7 = 1.0 and a = 0.25. These parameters 
are in agreement with that given by Figures 1 and 2. The error bars of Figure 3 are also 
the range of the 10 subsamples. Therefore, the HI mass density distribution in the scale 
range of ~ 1 to 15 h^^Mpc can also be described by the log-Poisson hierarchy. However, the 
parameters /? and 7 are different from that of density field p(x), which has /? ~ 0.47, and 
7 ~ 1 (Liu & Fang 2008). That is, the density field Phi(x) is weaker intermittent and less 
singular than that of p(x) . 
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Fig. 3. — Intermittent exponent ^{p) of simulation samples of the mass density field of neutral 
hydrogen. The solid line is given by equation (8) with parameters a — 0.25, /3 = 0.60 and 
7 = 0.75. The error bars are the variance of ^(p) over 10 samples, each of which contains of 
1000 one-dimensional samples. 



4. Log-Poisson hierarchy of Lya transmitted flux 

4.1. Samples 

4.1.1. Observed data 

For observed samples of Lya transmitted flux, we use the high resolution, high signal to 
noise ratio QSO Lya absorption spectra of Jamkhcdkar (2002), Jamkhcdkar et al. (2003). 
The power spectrum and intermittency of this data set have been extensively and deeply 
analyzed (Pando et al. 2002; Feng et al. 2002; Jamkhcdkar et al.2003). It is useful for 
testing the log-Poisson hierarchy. 

This observational data set consists of 28 Keck High Resolution (HIRES) QSO spectra 
(Kirkman & Tytler 1997). The QSO emission redshifts cover a redshift range from 2.19 to 
4.11. The resolution is about 8 km s~^. For each of the 28 QSOs, the data are given in form 
of pixels with wavelength, flux and noise. The noise accounts for the Poisson fluctuations in 
the photon count, the noise due to the background and the instrumentation. The continuum 
of each spectrum is given by IRAF CONTINUUM fltting. 

The data are divided into 12 redshift ranges from z — 1.6 + n x 0.20 to 1.6 + (n + 1) x 
0.20 where n = 0, . . . , 11. In our analysis below, we use only the data in the range z — 2.4 
to 2.6, which is the same as our simulation data. The scale range is taken to be 1.54 to 12.3 
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h ^ Mpc, which is also about the same as simulation samples. The mean flux in this redshift 
range is (F) ~ 0.75 (Jamkhedkar, 2002). 

We use the same methods of Jamkhedkar et al. (2003) to deal with noises, metal lines, 
proximity effect, and bad data chunks. On average the S/N ratio of the Keck spectra is 
high. Most of the regions with low S/N are saturated absorption regions. Although the 
percentage of pixels within these regions is not large, they may introduce large uncertainties 
in the analysis. We should reduce the uncertainty given by low S/N pixels. The method is 
as follows. First, we calculate the SFCs (scaling function coefficients) of both transmission 
flux field F{x) and noise field n{x) with eji = J F{x)(f)jj{x)dx and ^fi = J n{x)(f)jj{x)dx, 
where 4'j,i{x) is the scaling function of wavelet on scale j and at position We then identify 
unwanted mode (j, /) by using the threshold condition |ej^/e^| < /. This condition flags 
all modes with S/N less than f. The parameter / is taken to be 3. We also flag modes 
dominated by metal lines. In order to easily flag data gaps, we set the flux at the gaps to 
be zero and the error to one, and do the same thing for pixels with negative flux. Finally, 
we skip all the flagged modes when doing statistics. With this method, no rejoining and 
smoothing of the data are needed. 



To simulate Lya transmitted flux F, we use the 10,000 one-dimensional samples given in 
§3.1. For each 1-D sample, the transmitted flux of Lya, F[z) in redshift space is calculated 
with F{z) — exp[— t(2;)], where t{z) is the optical depth deflned as 



in which ctq is the effective cross-section of the resonant absorption; H is Hubble constant 
at the redshift of the sample, nui{x) is the number density of neutral hydrogen atoms. The 
Voigt function is V[z — x — v{x), b{x)] — eyiY){—[z — x — v{x)Y/h'^{x)}, and h{x) being 

the thermal broadening. To has a proper comparison between simulation and observation, we 
add noises rii — GxAiio each pixel i, where G is randomly sampled from standard normal 
distribution and Ai is the noise level of pixel i. We then take the same data reduction of 
noised modes as observed data. 



4.1-2. Simulation samples 




(20) 
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4.2. Redshift distortion and /^-hierarchy 

The velocity field v{x) and thermal broadening b{x) of eq.(20) will lead to the deviation 
of the statistical properties of In F{z) = —t{z) from n-m^x) or puiix). The field In F{z) is 
no longer to show all features of log-Poisson hierarchy. We should study which properties of 
the log-Poisson hierarchy can still be seen with Lya transmitted flux. 

Because the scale free range is larger than the Jeans length, the effect of thermal broad- 
ening would be small in this range. In this case, eq.(20) can be approximately as 



In F{z) — t{z) = -j^ J nui{x)S[z — x — v{x)]dx. 



(21) 



This relation actually is a mapping from a physical space fleld nm{x) to redshift space fleld 
t{z), which is the same as that used in the redshift distortion of galaxy distribution. It has 
been shown that, with the DWT variables, the redshift distortion of eq.(21) can be estimated 
by (Yang et al 2002) 

STr,l = mrSpr,l, (22) 

where the DWT variables Srr^i are given by 6Tr^i = f T{x)ipj^i{x)dx = J[— In F{x)]ipj^i{x)dx. 
The redshift distortion factor 5%. depends on the DWT power spectrum of the velocity field 
v{x) on scale r = . 



„" -0.1 



-0.2 




-0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 
ln[Fp(r)/F2(r)] 

Fig. 4. — The /9-hierarchy of observed sample of the Lya transmitted fiux at redshift z = 
2.4 — 2.6 and physical scale range from ~ 1 to 15 h~^ Mpc. The error bars are given by the 
maximum and minimum of bootstrap resampling. 

Because the average of (...) in the structure function eq.(l), or Sp{r — — {\Spj^i\P), 

is only over on modes I, the redshift distortion factor ^Rj. does not involve in this average. 
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Thus, the function Fp{r) of 5Tr,i will be different from the function Fp{r) of Spr,i by a factor 
yir- Thus, both Fp^i{r) / F^i^r) and Fp{r) / F2{r) of STr^i do not contain the redshift distortion 
factor yir- Thus, Fp+i{r) / F^i^r) and Fp{r) / F2{r) of STr^i should also satisfy the /9-hierarchy 
eq.(12) as that of field pni- That is, the /3-hierarchy is not affected by the redshift distortion. 

Figure 4 presents the /3-hierarchy of observed transmitted flux F{z). It is very well 
fitted by a straight line for data on the scale range from ~ 1 to 15 h~^ Mpc. It yields 
(3 = 0.67 ± 0.02. Figure 5 shows the /3-hierarchy of observed transmitted flux, which is the 
same as Fig. 4, and the noised simulation samples F. We see that the both real and simulation 
samples arc well coincident. The simulation samples yield /? = 0.66 ± 0.02. Therefore, both 
real and simulation samples are well /5-hierarchical, and the numbers of (3 are well consistent. 
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ln[Fp(r)/F2(r)] 
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Fig. 5. — The /^-hierarchy for a. observed data (square) as Fig. 4; b. simulation sample 
(circle) of the Lya transmitted flux at redshift z = 2.5 and Gaussian noises are added with 
the same level as real data. 

The data points of Figure 4 are scattered along the straight line Fp_^i[r) / F3[r)-Fp{r) / F2{r) , 
while the points of Figure 5 are clustered. This is due to all simulation samples have the 
same length, while for real samples consisting quasars with different redshift, we take only 
the sections of transmitted flux, which are in the redshift range z ~ 2.4 - 2.6. 



4.3. High order moments 

The statistics of eq.(14) is based on the ratio between high order moment {Sp'^.^) and 
{^PiY ) both of which have the same order of dp^. When p = 2, the ratio actually is kurtosis. 
The moment ratio of (5t^^) to {Sr^y, obviously, is independent of the redshift distortion 
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factor 9^. Therefore, the ratio of {5t^^) to {5t^Y has to satisfy the same property as {5pl^) 
to {SplY- Thus, for a given p, the relation of ln(((5rr)^^)/ {{SrrYY and Inr has to be a straight 
hne with slope — 7[p(l — — (1 — — 13) [eq.(15)]. Since the parameter j] is already 

determined by the /^-hierarchy, the test here is whether we can find one parameter 7 to fit 
the slops of {St^^) I {St'^Y vs. In r for different p. 



sr 

A 

CM 




r(h''Mpc) 



Fig. 6. — ln[((5rr)^^)/ ((^tv)^)^] vs. r of real data in the scale range ~ 1 - 15 h~^ Mpc, and p 
is taken to be 2 (bottom) and 3 (top). The solid lines are given by the least square fitting, 
which yield slopes 0.54 ± 0.19 (p = 2), and 1.3 ± 0.4 {p = 3). The error bars are given by 
the maximum and minimum of bootstrap resampling. 

The result is presented in Fig. 6. It shows the relation of ln{{ST^P) / {{STr^Y) vs. Inr 
for real data with p = 2 and 3, i.e. the statistics arc of the order of 4 and 6. The slopes 
of the fitted straight lines are 0.54 ± 0.19 for p = 2, and 1.3 ± 0.4 for p = 3. Thus, using 
p = 0.67, we have 7 = 0.59 ± 0.20 for p = 2 and 7 = 0.58 ± 0.22 for p = 3. In spite that the 
statistical errors are still large, one can already see that different p-lines yield about the same 
parameter 7. In other words, the moment statistics passed the test of log-Poisson hierarchy. 

The number of 7 given by the real data of the transmitted fiux is less than the number 
7 = 0.75 of Phi field (§3.2). Therefore, the field of transmitted fiux of real samples is less 
singular than the field pm- This is reasonable if considering the real data are noised. 



4.4. Log-Poisson hierarchy and log-normal model 

In §4.2 and 4.3., we have found that the log-Poisson parameters of the transmitted 
flux should he P — 0.67 ± 0.02 and 7 — 0.58 ± 0.22. This result can be used to com- 
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pare the log-Poisson hierarchy and log-normal model. We now consider the moment ratio 
ln[{ST^P) / {{5Tryy] as a function of p on a given scale r. 




Fig. 7. — ln{Sr^P) / {Sr^y as a function of p for real data of Lya transmitted flux at the 
range 2; ~ 2.4 - 2.6 on physical scales r = 1.54 (nabla), 3.07 (triangle), 6.14 (circle)and 12.3 
(square) h~^ Mpc. The curves are from the fitting of log-Poisson model. For clarify, wc show 
only the error bars for data points of r = 1.54 h~^ Mpc. The errors of other r are about the 
same level as r = 1.54 h~^ Mpc. 

Figure 7 presents the {6tI^'^) / {Sr^y as a function of p for real data on scales r = 1.54 
(nabla), 3.07 (triangle), 6.14 (circle)and 12.3 (square) h"^ Mpc. p is from 1 to 4. The 
p-dependent curves shown in Figure 7 are given by eqs.(14) and (15), in which the fitted 
parameters are f3 = 0.67 ± 0.02 and 7 = 0.55 ± 0.10 The const in eq.(14) is determined by 
\n[{5T^P)/{5T^)P] = when p = 1. 

Although the observed data points of Figure 7 actually are the same as Figure 6, we see 
that the observed data points show large scatter in Figure 6, but almost no scatter in Figure 
7. This is because Figure 6 gives moment ratio as a r function, while Figure 7 shows the 
p-dependence of the moment ratio. For a given r, the Gaussian noise will yield a moment 
ratio given by cq.(13), which is a smooth function of p. Therefore, a Gaussian noisy sample 
should not cause scatter with respect to p. On the other hand, Gaussian noise generally is 
not a smooth function of r. It yields the scatter of Fig. 6. 

Figure 7 shows that the p-dependcncc of moment ratio ln[(5r^^)/((5r^)P] is significantly 
dependent on scale r. Therefore, it can not be fitted by a Gaussian field, for which moment 
ratio is r-independent [eq.(13)]. The field of 6t is non-Gaussian. The smaller the scale r, 
the stronger the non-Gaussianity. 
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Fig. 8. — The same as Figure 7, but the curves are from the fitting of log-normal model. 



The log-normal model of the transmitted flux of QSO Lya absorption spectrum is 
very successful to explain various low-order statistical features of Lya forests (Bi 1993; Bi 
& Davidsen 1997). The log-normal model also predicts that the transmitted flux is non- 
Gaussian and intermittent. For the moment ratio, the log-normal model yields (Pando et al 
2002) 



(23) 



where (T^(r) is the power spectrum of the field. For each r, one can fit eq.(24) to observed 
points with o'^(r). The results arc plotted in Figure 8. It shows that, if we try to give a 
good fitting of eq.(24) with real data at orders p < 2, the p-curvcs of cq.(24) always give a 
large deviation from the real data at p > 2. This deviation cannot be reduced with selecting 
(7^(r). This is because the deviation is from the p^-dependence of {5t'^^^) / {5t'^Y when p is 
large. The p^-dependence cannot be reduced with the parameter (T^(r). 

On the other hand, at high p, the log-Poisson model gives {St^'') / {Sr^y oc p. The 
increasing with p is then consistent with observation. In Figs. 7 and 8 we show the error 
bars for data points of r = 1.54 h~^ Mpc. Although the errors are large at high p, the 
result is clearly consistent with p-dependence, and unfavor the p^-dependence. Therefore, 
the higher order statistics of the Lya transmitted flux is effective to discriminate between 
the log-Poisson and the log-normal model. The log-normal model yields too strong non- 
Gaussianity. This point actually has already been mentioned in the study of turbulence (e.g. 
Frisch 1995). 
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4.5. Scale-scale correlation 

The last statistical feature used to test the log-Poisson hierarchy is the scale-scale correc- 
tion. Similar to the statistics of high order moment (§4.3), the ratio of scale-scale correction, 
{St{!^St^^) / {5t^^){5t^^) is independent of the redshift distortion factor 9^^- The field r should 
show the scale-scale correlation as pni eqs.(r7)-(19). 



1 - 
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Fig. 9. — Scale-scale correlation of observed data of Lya transmitted flux at the range z ~ 2.4 
-2.6 for p — 2 and r2/ri = 2. The error bars are given by the maximum and minimum of 
bootstrap resampling. 

Before doing this test, we would emphasize that the scale-scale correlation is statistically 
independent of the statistics of high order moments. For instance, one can construct a field, 
which shows Gaussian distribution in terms of its one point distribution of Spj^i, but highly 
scale-scale correlated between variables Spjj with different j (Pando et al 1998, Feng & Fang 
2000). The clustering of cosmic large scale structure in the nonlinear regime essentially is due 
to the interaction between Fourier modes on different scales. Therefore, cosmic clustering 
will definitely leads to scale-scale correlation. Scale-scale correlation is also very effective to 
distinguish various hierarchy cascade models (Pando et al. 1998). 

Figure 9 presents the p = 2 scale-scale correlation of observed data, in which the ratio 
r2/ri is fixed to be equal to 2. The slop of the fitting straight line in Fig. 9 should be 
given by eq.(19). Since eq.(19) depends only on the parameters /? and 7, both of which have 
already been determined in §4.2 and 4.3. There is no free parameters in the fitting of Fig. 9. 
If we doing a straight line, we found the parameters /3 = 0.67 and 7 = 0.43±0.12. The value 
of P is the same as that of §4.2, while the value of 7 is smaller than that of §4.3, but the 
deviation is not larger than 1-a. The scale-scale correlation is more sensitive to the quality 
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of the data, as it is the correlation between modes on different scales. The result of Figure 9 
is basically consistent with the scale-scale correlation predicted by the log-Poisson hierarchy 

5. Discussion and conclusion 

Nonlinear evolution of mass and velocity fields is a central problem of large scale struc- 
ture of the universe. The clustering of the cosmic baryon fluid, governed by the Navier-Stokes 
equation in gravitational field of an expanding universe, has to be self similarly hierarchi- 
cal in the scale free range in which the dynamical equations and initial perturbations are 
scale- invariant. 

The log-Poisson hierarchical clustering sketches the nonlinear evolution of cosmic baryon 
fluid in the scale free range. If the initial density perturbations are Gaussian, and their power 
spectrum is given by power law P{k) oc A;", the structure functions initially have to be Sp{r) oc 
In the regime of linear evolution, the structure functions will keep to be Sp{r) oc r^^^p\ 
and the intermittent exponent is ^i{p) = —pa/2. According to the log-Poisson hierarchy 
scenario, the nonlinear evohition leads to the hierarchical transfer of the power of clustering 
from large scales to small scales. The structure function will become Sp{r) oc r^'^p'i+^niip) ^ 
where the nonlinear term of the intermittent exponent is C,ni{p) = — 7[p— (1 — — in 
which the parameters (3 and 7 are dimensionless. /3 measures the intermittency of the field, 
and 7 measures the singularity of the clustering. For Gaussian field, we have (3 — 1, and 
therefore, Cni{p) — for all order p. Since the onset of nonlinear evolution, the parameter P 
will gradually decrease, and the field becomes intermittent. 

With simulation samples, we found that the parameter P is decreasing with the decrease 
of redshift z. It means that the field is weakly intermittent at earlier time, but strong 
intermittent at later time (Liu & Fang 2008). Although ^ni{p) 7^ 0, the nonlinear evolution 
keeps the cosmic baryon fluid to be scale-invariant. We showed that the mass density field 
of neutral hydrogen fiuid in the scale free range is also well described by the log-Poisson 
hierarchy in spite of the neutral hydrogen fraction of the baryon fiuid is not constant in 
space. This is because the UV ionization photon is assumed to be uniform, and it does not 
violate the scale invariance of this system. However, the number of (3 of neutral hydrogen 
is found to be less than that of total baryon fluid. Therefore, the neutral hydrogen is less 
intermitted. 

The Lya transmitted flux of quasars Lya absorption spectrum is considered to be ef- 
fective to probe the mass and velocity flelds of neutral hydrogen. However, the redshift 
distortion of the velocity held leads to the deviation of the held of the Lya transmitted flux 
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from the neutral hydrogen field. The transmitted flux does not satisfy all predictions of 
log-Poisson hierarchy. Fortunately, the effect of radshift distortion is approximately negligi- 
ble for some log-Poisson hierarchical predicted features. Thus, we can test the log-Poisson 
hierarchy with quasars Lya absorption spectrum. 

Using high resolution and high S/N data of quasars Lya absorption spectrum, we show 
that all the non-Gaussian features predicted by the log-Poisson hierarchy, including the 
hierarchical relation, the intermittent exponent, the ratios of different moments, and the 
scale-scale correlation, are consistent with observed samples. The observed samples of the 
transmitted flux yield the same intcrmittence parameter /? as that of neutral hydrogen field 
produced with hydrodynamic simulation of the concordance ACDM universe. Our result 
shows that the intermittent exponent i{p), or parameters /3 and 7, is effective to discriminate 
among models of nonlinear evolution. 

The log-normal model can well fit observed data on lower order statistics, but not good 
on higher orders. On the other hand, the log-Poisson model gives good fitting on lower order 
as well as higher order statistics. Therefore, a comparison between the log-Poisson model 
and log-normal model on lower order statistics will be able to find the relation between 
parameters of the log-Poisson and log-normal models. This relations would be useful to 
explain the parameters of log-Poisson model with well-known parameters in cosmology, as 
the parameters of log-normal model generally are known in cosmology. 

Recent studies have shown that the turbulence behavior of baryon gas can be detected 
by the Doppler- broadened spectral fines (Sunyaev et al. 2003; Lazarian & Pogosyan 2006). 
Although these works focus on the turbulence of baryon gas in clusters, the result is still 
applicable, at least, for the warm-hot intergalactic medium (WHIM), which is shown to 
follow the evolution of Burger's fluid on large scales (He et al. 2004, 2005). The last but 
not least, the polarization of CMB is dependent on the density of electrons, and therefore, 
the map of CMB polarization would provide a direct test on the non-Gaussian features of 
ionized gas when the data on small scales become available. 

We thank Ji-Ren Liu for his contributions in the early stage of this project. Y. Lu is 
supported by China Scholarship Council. This work was partially supported by US NSF 
AST 05-07340. 
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